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ABSTRACT 

O 

^ . The perturbation method is apphed to numerical solution of the Lane-Emden 

r\ • Equation (LEE)of arbitrary index n, and the global parameters of polytropes are 

^ . found as function of polytropic index n. 

in 

>• ' 1. Introduction 

o 

m 

The basic differential equation of internal structure of stars, the Lane-Emden equation 

^ . (hereafter LEE) of index n, is solved analytically only for three values of = 0, 1, 5. For other 

'nI" . values of n, LEE can be solved only numerically. Recently, using the perturbation method, 

. SK ( Seidov & Kuzakhmedov, 1978 (SK78)) had presented the new analytical solutions of the 

P_|. LEE for index n only slightly differing from 0, 1, and 5, see also Seidov (1978a, b, 1979a,b, 

Q : 2004); Jabbar (1984); Caimmi (1987); Horedt (1987, 1990); Medvedev & Rybicki (2001). 

+i . In this paper I present the numerical perturbation method for solving the LEE, see also 

^ : Seidov (2004). 
> 

X 

■ 2. Basic equation 

The basic equation is LEE of index n : 

LE = ^^(x'4-]-^ LE[y] = -y-- y(0) = 1, 2/'(0) = 0, 



where we introduced the Lane-Emden differental operator LE. 

We look for solution y{x) in the interval [0,X] such that y{X) = 0. 

Three classical analytical solutions of equation (1) are (see e.g. (Chandrasekhar 1957)): 

n = 0, y = l-^x^, X = Vg, /i = 2v/6, pc/pm = l; (2) 



n = l, y = X = TT, /i = TT, pc/pm = ttVS; (3) 

X 
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n = 5, y = (1 + - x^) 1/2^ X ^ oo, 11 = VS, pc/pm oo. (4) 

In these equations, p = —X'^y'{X), pd Pm = -^^/S/^; X, p are dimensionless radius and 
mass, and pd pm is the central-to-mean density ratio. 



3. The perturbation method 

Consider equation (1) as ODE depending on parameter n, then assuming n = uq + 5 
with 5 as a small parameter, 5 -C no, (or 5 ^ 1, if no = 0) we expand the r.s. of equation 
(1) to the second order of 5 (for the sake of brevity we omit index at n): 

y^yo + Syi + S^y2; yi;+^ = ho + 5 ho + 

ho = Vo; hi = ny'^'^ yi + y"^ Int/o; 

1 1 
/i2 = - n (n - 1) y'^-^ yj + y'^'^ {yi+ny2+n yi In yo) + -y^ In' yo- (5) 

Prom equations (1,5) we have three coupled ODEs for three functions yo, yi, y2: 

LE[yo] = -ho, yo(0) = l, y;(0) = 0; (6) 

LE[yi] = -hi, yi(0)=y;(0) = 0; (7) 

LE[y,] = -h,, y2(0) = y;(0) = 0. (8) 

Initial conditions in equations (6,7,8) are defined by the form of series expansion of the 
solution of LEE of arbitrary n at x = (see further, formulas (16 - 22): 

y= l-^x2 + ^x^ + .... (9) 

Writing n = no + S, expanding equation (9) to the second order of 5, we have the series 
expansions for functions yi, y2 at x = 0: 

1 4 5-16n 6 /.n^ 

8 (5 183 (-1 + 2 n) o 
= 9.9! 

Note that the series expansion for yi and y2 were given in Seidov (2004) only for case n = 0. 
Before solving eqs (6,7,8), I'd like to mention that the vahdity of the approach used is 
discussed briefly in SK78, and in Seidov (2004). For n = 0, 1, and 5 the ana/y^ica/ solutions 



-3- 



for functions yi were presented in SK78; also the second approximation at the case of n = 
is partly given in Seidov (2004). Here we solve numerically equations (6-8). We remind 
that I/O is the "basic" solution of LEE and we refer to yi and y2 as the "perturbed" LEE 
solutions of the first and second order. I widely used MATEMATICA^s (Wolfram 1999) 
function NDSolve with suitable options. 

4. System of ODEs 

For a better numerical accuracy, we introduce three additional functions, zq, zi, and Z2: 
zo = -X yo, zi = -X y^, Z2 = -x y^, (12) 
and rewrite the three equations (6-8) as a system of six differential equations: 

y'o = -^o/a;^ 4 = ^0, (13) 

y'l = -zi/x'^, z[ = -x^ hi. (14) 

y2^~Z2/x^, Z2^—X^h2. (15) 



5. Series solution of zero-th order 

Again, for a larger numerical accuracy, we use the series solution of LEE. In Seidov & 
Kuzakhmedov (1977), the method of accurate series solution of LEE is given. Using the 
formulas of Seidov & Kuzakhmedov (1977), we present solution of LEE (13) at a; = in the 
form 

i=12 i=12 

yo = l + 5] -2o = -2 5^ ia^x^^+S (16) 

i=l 1=1 

with coefficients Oj as follows: 

_ 1 _n _n{5-8n) _ n (70 - 183n + 122n2) 

"^-"6' "'-5!' 3-7! ' 9T9! ' 

_ n(3150 - 10805 n + 12642 - 5032 n^) 
45T11! ' 
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_ n (138600 - 574850 n + 915935 - 663166 + 183616 n^) 

135^13! ' 
ay = n (21021000 - 101038350 n + 199037015 

-200573786 n-'^ + 103178392 - 21625216 n^)/945 ■ 15!; 

as = n (1891890000 - 10267435500 n + 23780949500 - 30057075285 

+21827357636 - 8618115372n^ + 1442431856 n^)/2835 ■ 17!; 

ag = n (675404730000 - 4066235428500 n+ 10740122081500 

-16120795594195 + 14830640277988 - 8348507232868 

+2657923739344 - 368552598784 n^)/25515 • 19!; 

aio = n (171102531600000- 1128186384570000 n 

+3329284073314500 - 5740042719521900 

+6317195348852735 - 45381 14873629364 

+207492591889115671^ - 551199819782480 

+65035924972928 n^)/127575 • 21!; 

an = n (118574054398800000 - 8479530565991 10000 n 

+2754994980587692500 - 5335484162711174500 

+6782008348777403475 - 5860922969087284308 

+3438918097715059380 - 1319254687791147504 

+299840088682556928 - 30720693974199296 n^)/ 1403325 • 23!; 

ai2 = n (27272032511724000000- 209899877314257900000 n+ 742585473289204545000 

-1589853990586539282500 + 2279636465710370388750 - 2285217511971127632065 

+1620103707989338077938 - 801095938682391176900 + 264081052577164986584 

-52342890902954850528 + 4731477379473053696 n^°)/4209975 • 25!.(17) 

We remind that the general recurrence relation for coefficient a^ is (Seidov & Kuzakhmedov 
(1977)): 

^ i=m 

Qm+i = —tjt: r^yZ(^'^ + ^ ~ m)(m + 1 -i)[3 + 2(m - i)]aj a^+i_j. (18) 

m[m + lj(2m + 6) ^—^ 
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6. Series solution of the first order 

Writing n = no + 5, expanding equation (17) to the first order of 6, we have, from (16), 
the series expansions for solutions yi, z\ of the equations (14) at a; = 0: 

i=12 j=12 

7/1 = 1 + ^i = -2 5^^6,a;2^+\ (19) 

1=2 1=2 

with coefficients hi as follows: 

62 = 1/5!; 63 = (5- 16n)/3-7!; 64 = (70 - 366 n + 366 n^)/9 ■ 9! 

65 = (3150 - 21610 n + 37926 - 20128 n^) /45 ■ 11! 
fee = (138600 - 1149700 n + 2747805 r? - 2652664 + 918080 n^)/135 ■ 13! 

67 = (21021000 - 202076700 n + 597111045 - 802295144 
+515891960 - 129751296 n^)/945 • 15!; 
68 = 4(472972500 - 5133717750 n+ 17835712125 - 30057075285 
+27284197045 - 12927173058 + 2524255748 n^)/2835 • 17!; 
69 = 4(168851182500 - 2033117714250 n + 8055091561125 
-16120795594195 + 18538300347485 
-12522760849302 n^ + 4651366543852 - 737105197568 n^)/25515 • 19!; 
610 = (171102531600000 - 2256372769140000 n + 9987852219943500 

-22960170878087600 + 31585976744263675 
-27228689241776184 + 14524481432238092 
-4409598558259840 + 585323324756352 n^)/127575 • 21!; 
611 = (118574054398800000- 1695906113198220000 n 
+8264984941763077500 - 21341936650844698000 
+33910041743887017375 - 35165537814523705848 
+24072426684005415660^6 - 10554037502329180032 
+2698560798143012352 - 307206939741992960 n^)/1403325 ■ 23!; 

612 = 2 (13636016255862000000- 209899877314257900000 n 
+ 1113878209933806817500^2 - 3179707981173078565000 
+5699091164275925971875 - 6855652535913382896195 
+5670362977962683272783 - 3204383754729564707600 
+1188364736597242439628 - 261714454514774252640 n*^ 

+26023125587101795328 n^°)/ 4209975- 25!. (20) 



-6- 



7. Series solution of the second order 

Again, as in section 6, writing n = uq + S, expanding equation (17) to the second order 
of 5, we have get the series expansions for solutions y2, Z2 of the equations (15) at a; = 0: 

i=12 i=12 

y2 = l + ^Qx2\ ^2 = -2^ZQx2^+\ (21) 

i=3 1=3 

with coefficients Cj as follows: 

C3 = -l/3-7!; C4 = 61(-l + 2n)/3-9!; 
C5 = (-10805 + 37926 n - 30192 71^)745 • 11!; 
C6 = (-574850 + 2747805 n - 3978996 + 1836160 n=^)/135 • 13!; 

C7 = (-101038350 + 59711104571- 1203442716 
+1031783920 - 324378240 n^)/945 • 15!; 
Cg = 2 (-5133717750 + 35671424250 n - 90171225855 
+ 109136788180 - 64635865290 n^+ 15145534488 n^))/2835 ■ 17!; 
Cq = 2 (-2033117714250 + 16110183122250n- 48362386782585n^ 

+74153201389940 n-^ - 62613804246510 n'^ 
+27908199263112 - 5159736382976 n^)/25515 ■ 19!; 
cio = 2 (-564093192285000 + 4993926109971750 n 
-17220128158565700 + 31585976744263675 
-34035861552220230 + 21786722148357138 
-7716797476954720n^ + 1170646649512704 n^)/127575 • 21!; 
cii = 2(-141325509433185000+ 1377497490293846250 n 

-5335484162711174500 + 11303347247962339125 
-14652307422718210770 + 12036213342002707830 
-6156521876358688352 + 1799040532095341568 rJ 
-230405204806494720 n^)/467775 • 23!; 
= (-209899877314257900000 + 222775641986761363500071 
-9539123943519235695000 71^ + 2279636465710370388750071^ 
-34278262679566914480975 7i^ + 34022177867776099636698 
-22430686283106952953200 + 9506917892777939517024 ri^ 
-2355430090632968273760 + 260231255871017953280 n^) / 4209975 • 25!. (22) 
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8. Cases with analytical solutions 

We first consider three cases, for wliicli tlie zero-tli and first approximations liave exact 
analytical solutions. 



8.1. n=5 case 



We first mention that in this case coefficients Oj can be found from the very simple 
recurrence relation: 

n = 5: arn = ^ (7: - m] Grn-i, m>l, ao = 1, (23) 
3m \2 / 

which can be used e.g. for checking formulas (16, 17). Then, we calculate, for n — 5, values 
of yo and zo according to (16) and (17) at point Xi = 1/100. With using MATHEMATICA 
it can be done exactly (without any numerical error). Then we compare these values with 
exact analj^ical values of yo and Zq at Xi — 1/100 (see (4)): 

(l + x,V3)-V2 

~0 l.j.Ki 

Hence we have some 60 digits of accuracy in the initial values for LEE (13). Now we 
solve the system (13-15) numerically using NDSolve procedure of MATHEMATICA, from 
X — Xi — 1/100 till X — Xf — 10 with options: Method — > RungeKutta, AccuracyGoal — > 
Infinity, PrecisionGoal — > 32, WorkingPrecision — > 24, MaxSteps 50000, and compare 



- 1 = 9.72065376 • lO"'*", 

- 1 = -2.5273216 • 10"™. (24) 
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the numerical values with values from analytical expressions when available 



n — 5 Xf — 10 





, NDSolve 


0.17066 40371 96572 26797 786 


yo{xf] 


, Analytic 


0.17066 40371 96572 28860 143 




, NDSolve 


1.65693 23999 66721 38805 555 




, Analytic 


1.65693 23999 66721 24855 754 




, NDSolve 


0.096926 10440 49489 40514 448 




, Analytic 


0.096926 10440 49489 36327 188 


Zi{Xf] 


, NDSolve 


-0.12538 65905 38367 83413 8922 




, Analytic 


-0.12538 65905 38367 81528 4516 


y2{xf] 


, NDSolve 


-0.0131 49071 88483 68785 




, NDSolve 


0.0001 04233 84608 82482. 



Final results arc of 15-16 digits accuracy. Here "analytical" values for yl and zl are given 
according to analytical solution SK78: 

n — 5 : V — arctan(-i^), 
yi = — — . . , ( sin(2i/) — - sin(4i/) + 31/ cos(4i/) — 3 [2 sin(2i/) + sin(4i/)] ln[cos(i/)] j , 

4to SlHl 1 \ 4t J 



^dyi 



z, = -x'^. (26) 
ax 



8.2. n=l case 

In the case of n = 1 we first mention that coefficients are of very simple form: 

n=l: «.= (^-^, (27) 

which can be used e.g. for checking formulas (16, 17). Then, we calculate, for n — 1, values 
of yo s-nd Zq according to (16) and (17) at point Xi — 1/100, and compare these values with 
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exact analytical values of yo and zq at Xi — 1/100: 



y^^yi^;'-' - 1 = 9.183841796 • lO'^S 
sm{xi)/Xi 



^o(3^i) 16,17 



sin(xi) — Xi cos(xi) 



- 1 = -2.3877590 • 10 



-81 



(28) 



Now we have some 80 digits of accuracy in the initial values for LEE (13). Then we solve 
the system (13 - 15) numerically using NDSolve procedure of MATHEMATICA, from x = 
Xi — 1/100 till X — Xf — TT with the same options as for case n — 5, section 8.1, and compare 
the numerical values with values from analytical expressions when available 



n 



X = 7: 



I/O (X), NDSolve 
I/O (X), Analytic 



^o(X), NDSolve 
2;o(X), Analytic, tt 



1/1 (X), NDSolve 
l/i(X), Analytic 
^Si(27r) + iln(27r)-| 

^i(X), NDSolve 
Zi (X) , Analytic 



-6.552265490 ■ lO'^^ 


3.14159 26535 89793 22713 584 
3.14159 26535 89793 23846 264 

0.28179 14499 02207 82015 93981 

0.28179 14499 02207 82012 97922 

-1.02925 45490 49509 89738 663 



|Si(27r) + f [Ci(27r) + ln(27r) - 3 - EulerGamma] -1.02925 45490 49509 87940 655 



1/2 (X), NDSolve 
1/2 (X), Analytic 



^2 (X), NDSolve 
2:2 (X), Analytic 



-0.09455 03188 95873420 



-0.811718053186985 23026 



(29) 



Again, as n = 5 case, final results are of 15-16 digits accuracy, and though there is no 
guarantee that the similar accuracy will be in other numerical results, still we will use the 
same set of options for solving other cases of n. 
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8.3. n=0 case 

In the case of n = we first mention that coefficients are of very simple form: 

n = 0: oi = aj>i = 0; (30) 
while coefficients bi can be found from the series expansion of the known function yi (SK78): 

n = 0: = ^-4 + (2 + yo) ln(yo)+4^arctanh(^), (31) 

with (see equation (2)): 

6 

and 

Z/i(V6)=41n2--. (32) 

Then we solve, for n = 0, the system (13-15) numerically using NDSolve procedure of 
MATHEMATICA, from x = = 1/100 till x ^ X ^ with the same options as for cases 
n — 5, 1, and compare the numerical values with values from analytical expressions when 
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available: 



n= 

yo(^),NDSolve 
yo(^), Analytic 



Zo(X),NDSolve 
;2o(X), Analytic = 2^6 

l/i(X),NDSolve 
7/1 (X), Analytic = 4 ln2 - | 

Zi(X),NDSolve 
Analytic = 1^/6 (3 In 2 - 4) 

l/2(X),NDSolve 
1/2 (^), Analytic) = 
: i(413-2l7r2-402 In2 + 144 1n2 2) 

^2(A:),NDSolve 
2^2 (^), Analytic 



X = y/6 

1.530418314373 • 10"^^ 


4.89897 94855 66356 19639 45681 49410 81340 
4.89897 94855 66356 19639 45681 49411 78278 



0.43925 53889 0644790287 
0.43925 53889 0644790433 

-6.272517658760946466 
-6.272517658760954 355 

-0.4135138937125 

-0.4135138893 059 

-2.59850 • 10^2 
oo 



(33) 

We see that numerical accuracy is drastically decreasing for the larger degrees of pertur- 
bation: while yo{X), and Zo{X) are calculated with some 30-31 correct digits, yi{X), and 
^i(X) are calculated with 15-16 correct digits, and y2{X), and Z2{X) only with 8-9 correct 
digits. 



9. Non-analytical cases 

For other values of n, there is no analytic solution even for yo and there is no possibility 
to check numerical results. 

Moreover we should calculate numerically LEE to find accurate values of the first zero X. 
Also we give comparison of our numerical values of X and zo{X) with two most accurate 
results known in literature (Jabbar (1984),Horedt (1990)). 



9.1. n=l/2 case 



n = 1/2 



X, NDSolve 2.75269 80540 64991 

X, Jabbar (1984) 2.75269 80541 

X, Horedt (1990) 2.75269 8 



yo{X), NDSolve -1.57415358 • 10" 

1/0 (X), Analytic 



Zo{X), NDSolve 3.78865 11848 84005 74354 209 

Zo{X), Horedt (1990) 3.78865 11652 97791 



yi{X), NDSolve 0.34123 20401 31647 34451 



zi (X), NDSolve -1.62970 77128 59769 520 



2/2 (X), NDSolve -0.1468213064 84602 959 



Z2(X), NDSolve 



-3.7185797123 08979172 



9.2. n=3/2 case 



n = 3/2 




X, NDSolve 


3.65375 3736219119 


X, Jabbar (1984) 


3.65375 37362 


X, Horedt (1990) 


3.65375 4 


7/0 (X), NDSolve 


6.8926652958 • 10"^^ 


yo(^), Analytic 





Zo(X), NDSolve 


2.71405 5120108646 


^o(^), Horedt (1990) 


2.71405 57437 


yi(X), NDSolve 


0.240670919140952 32 


^i(X), NDSolve 


-0.71009 78273 59729 


t/2(^), NDSolve 


-0.0665618185 10763 63 


Z2(^), NDSolve 


0.2418612563 9746212 



9.3. n—2 case 

In this case we first mention that recurrecy relation for coefficients a^, may be obtained 
in the more simple form than in (17). This relation, according to Seidov (1979a), is: 

^ k=i 

"-+^^- (2z + 2)(2z + 3) ^"^^^^ ^^^^ 

Prom equation (36) it follows that the series is sign-alternating, the result not evident from 
(17). In general, it can be shown, that coefficients of serial solution for LEE of integer 
index n is sign-alternating. 
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We present results of numerical solution of the system (13-15). 



n = 2 




X ND Solve 


4 352874595946124 


X,Jabbar (1984) 


4.3528745959 




4 '^'i987 




SQ1 SQ9fi'i1 074 • 1 n~^'^ 


yo(^), Analytic 





^o(^),NDSolve 


2.41104 60120 96894 


;2o(X),Horedt (1990) 


2.41104 73856 


yi(X),NDSolve 


0.20990106855590432 


^i(X),NDSolve 


-0.51626996767525199 


7/2(^),NDSolve 


-0.049552916011399585 



Z2{X), NDSolve 0.15512703899482025 
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9.4. n=5/2 case 

We present results of numerical solution of the system (13-15). 



n = 5/2 



X, NDSolve 


5.355275459010824 


X, Jabbar (1984) 


5.3552754459 


X, Horedt (1990) 


5.355275 


7/0 (^), NDSolve 


-3.385331 • 10-1^ 


yo(^), Analytic 





zo(^), NDSolve 


2.18719 9565517079 


^o(^), Horedt (1990) 


2.18719 90907 


2/1 (X), NDSolve 


0.18557608273853998 


^i(X), NDSolve 


-0.38690411187876034 


y2(^), NDSolve 


-0.03834028384550369 



Z2{X), NDSolve 0.10810493285592643 
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9.5. n=3 case 



n — 3 




X, NDSolve 


6.89684 86193 7696 


X, Jabbar (1984) 


6.89684 86194 


X, Horedt (1990) 


6.89684 9 


I/O (X), NDSolve 


1.7861333071518752-10-^^ 


Vn(X), Analytic 





zo(^), NDSolve 


2.01823 59509 66228 7 


Zo(X), Horedt (1990) 


2.01823 62876 


yi(X), NDSolve 


0.1654768294 39454 83 


^i(X), NDSolve 


-0.29335 99255 39480 26 


7/2 (^), NDSolve 


-0.03053 51512 78657177 


Z2(X), NDSolve 


0.08152 7581187124 82 
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9.6. n=7/2 case 



n = 7/2 




X, NDSolve 


9.53580 53442 4479 


X, Jabbar (1984) 


9.53580 53443 


X, Horedt (1990) 


9.535805 


I/O (X), NDSolve 


1.2322310495888572 • 10"^^ 


Vn(X), Analytic 





zo(^), NDSolve 


1.89055 70934431164 


Zo(X), Horedt (1990) 


1.89055 65987 


yi(X), NDSolve 


0.1481615824132385 


^i(X), NDSolve 


-0.21964831524439676 


7/2 (^), NDSolve 


-0.024934678004649848 


Z2(X), NDSolve 


0.06771308534193725 



9.7. n=4 case 



n = 4 




X, NDSolve 


14.97154 63488 38093 


X,Jabbar (1984) 


14.97154 63496 


X, Horedt (1990) 


14.97155 


7/0 (X), NDSolve 


9.088953198578 • 10"^^ 


yo(^), Analytic 





Zo(X), NDSolve 


1.79722 99144 3925 


^o(^), Horedt (1990) 


1.79723 08344 


yo(^), NDSolve 


0.0010711089 71589 365 


yi(X), NDSolve 


0.13250 00877 72273 87 


^i(X), NDSolve 


-0.15399 30629 85788 74 


y2(^), NDSolve 


-0.0209711409 23665 716 


^2 (^), NDSolve 


0.06552 46752 2523159 
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9.8. n=9/2 case 



n = 9/2 




X, NDSolve 


31.83646 32446 95228 


X, Jabbar (1984) 


31.83646 32485 


X, Horedt (1990) 


31.83646 


l/o(X), NDSolve 


-1.6287136623021175 • 10" 


7/0 (X), Analytic 





zo(^), NDSolve 


1.73779 88676 66032 3 


Zo(X), Horedt (1990) 


1.73779 86022 


yi(X), NDSolve 


0.11719722297815073 


^i(X), NDSolve 


-0.07993661318489896 


1/2 (X), NDSolve 


-0.018685417777349604 


Z2(X), NDSolve 


0.08780938436898907 



10. Global parameters 
10.1. Radius 

Radius of the polytrope is proportional to X, the first zero of y{x). In our second order 
approximation in S, we write: 

X^Xo + SX. + S-'X^, yiX) = 0, y{X) = y^iX) + Sy.iX) + y^iX), (43) 

from where we get: 

^^^MXol, _ y2(Xo) + X, y[{X,) + 1/2 XI (Xq) _ ^^^^ 

^ -yoW ^ -yo(^o) 
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Let me remind the meaning of Xq, Xi, and X2: for each n, Xq is the dimensional radius, 
while Xi, and X2 define the derivatives of X{n) by n: Xq — X{n), Xi — dX{n)/ dn, and 
2X2 = d'^X{n)/dn^. 



10.2. Mass 

Mass of the polytrope is proportional to ji — —X'^y'(X). In our second order approxi- 
mation in 5, we write: 

H ^ IJ,o + 5 iJ,i + 5^ IJ,2, (45) 

from where we get: 

/io = -X^ y'oiXo), /ii = -X^ [Xi y^iXo) + y[{Xo)] - 2X, X, y^{Xo). (46) 

112 = -XI [X2 yo (Xo) + I XI yo'(Xo) + X^ yl{X^) + y^{X^)\ 
-2X^X^{X^yl{X^)+y[{X^))-y^{X^) {XI + 2X^X2). ^ ' 

Let me remind the meaning of /xq, A*i, and ^2- for each n, //q is the dimensional mass, 
while //i, and //2 define the derivatives of //(n) by n: //q = /^(^t-), A*i = d^{n)/ dn, and 
2^2 — d'^ ^{n)/ d-n? . 



10.2.1. n^O 



n = 0, 
^0 = 

y'^^{Xo) = -0.81649 65809 277260, 
y'^iXo) = -0.33333 33333 33333, 

y'oiXo) = 0, 

y^{Xo) = 0.43925 53889 064479, 

y[{Xo) = 1.0454196097934910777, 

y'^{Xo) = 70.1010152396 6198, 

y^(Xo) = -0.413513893712500 52051096, 

y'^{Xo) = 4.33084 54146 19353 • 10^\ 

Xi = 0.537975784794289 96 7933774, 

X2 = 0.12328 30846 88831 89 0159441, 

X{n = .1) = 2.50452 0152109495, (*2.5045449630 Jabbar*) 

X{n = -.1) = 

/io = 4.8989794855 66356196, 
= -3.04466 29499 952067, 
//2 = -2.59850 72489 991943 • 10^^ 
fx{n = .1) = 4.59451 31905, (*Horedt*) 
H{n = -.1) = -. 
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10.2.2. n^l/2 

n = 1/2, 

Xq = 2.75269805406499, 

yg(Xo) = -0.49999708294422937513122874, 

yl{Xo) = 0.3632777637699180575486, 

I/q"(Xo) = 3.7771115383 ■ 10^ 
1/1 (Xo) = 0.34123204013164497867521625893, 
y[{Xo) = 0.21507634824889584579L 
y[{XQ) = -2.723045849967156701 ■ 10^ 
1/2 (Xo) = -0.146821475444243995570334, 
l/2(Xo) = 929115.93812128759, 

Xi = 0.6824680618580861779508, (49) 
X2 = 0.169124728807637034976, 

X{n = .6) = 2.822636107538865, (*2.8226750739, Jabbar*) 
xln = .4) = 2.686142495167248, (*2.6861053263, Jabbar*) 
/io = 3.7886511848840057435396, 
/ii = -1.6297077071514047746998, 
1^2 = 0, 

= .6) = 3.6256804141688654, (only 1st appr!) 
(* 0.4560739 * 2.822675^ = 3.63376 61327, Horedt 
pi(n = .4) = 3.9516219555991463, only 1st appr! 
(* 0.5489336 * 2.686105^ = 3.96064 37923, Horedt *). 



10.2.3. n^l 



Xq = TT, 

y;,(Xo) = -0.31830988618379067039, 

yo(^o) = 0.202642367284675542812372, 

^^"(Xo) = 0.12480067956938'3.2195, 

yi(Xo) = 0.2817914499022078201593, 

y[{Xo) = 0.104285289178956941658, 

yl{Xo) = -0.348181526960602809649062, 

y2{Xo) = -0.09455031889587342081743776, 

y'^(Xo) = 0.0822442339327594290586, 

Xi = 0.88527394885719235129261239, 

X2 = 0.2424591740663775292855118679, 

X{n = 1.1) = 3.232544640216176, (* 3.2326084072, Jabbar*) 

X{n = 0.9) = 3.0554898504447374, (* 3.0554293447, Jabbar*) 

/io = 3.1415926535897932271, 

/ii = -1.02925454904951, 

fi2 = 0.4193306804206245, 

^(n = 1.1) = 3.0428605054890485, 

(*0.2911738 * 3.232608^ = 3.042694721493137, Horedt*) 

fx{n = .9) = 3.2487114152989505, 

(*0.3055429 * 3.0554292 = 3.24889082611291, Horedt*). 



-24- 



10.24- n^3/2 



n = 3/2, 

Xo = 3.653753736219119, 
y'Q{Xo) = -0.203301282638546737276897, 
y'^{Xo) = 0.111283516797123609328, 
y™(Xo) = -0.09137193541041363, 
yi{Xo) = 0.24067091914095230177155721, 
y[{Xo) = 0.0531911817234202858, 
y'^(Xo) = -0.02911591910998444970746, 
y2{Xo) = -0.066561818510763637968669, 
^^(Xo) = -0.018117061781318394264, 
Xi = 1.183814071497255432331730, 
X2 = 0.3658800713499040420899, 

X{n = 1.6) = 3.7757939440823436, (* 3.77590 47640, Jabbar*) 
X{n = 1.4) = 3.5390311297828925, (*3. (* 3.53892 66160, Jabbar*) 
/io = 2.7140551201086457355, 
/xi = -0.71009782735972906926, 
/i2 = 0.24186141001555074558, 
,_i(n = 1.6) = 2.64546 39514 728284, 

(* 0.185544* 3.7759052 ^ 2.64538 58927, Horcdt*)/i(n = 1.4) = 2.78748 3516944774, 
(*0.2225779 * 3.538927^ = 2.78756 657920, Horedt*). 

(51) 



10.2.5. n=2 



n = 2, 

Xo = 4.352874595946124, 
y'^{XQ) = -0.12724865113117523, 
yl{Xo) = 0.05846649074139796, 

^^"(Xo) = -0.040295089683388775, 
yi{Xo) = 0.20990106855590432, 
y[{Xo) = 0.027247367605845233, 

yl{Xo) = -0.012519252280422236, 2/2(Xo) = -0.049552916011399585, 
l/;(Xo) = -0.008187196082955055, 
Xi = 1.6495347234724416, 
X2 = 0.5888879298143976, 

X{n = 2.1) = 4.523716947591513, (* 4.5239262993, Jabbar*) 
X{n = 1.9) = 4.193810002897024, (* 4.1936146217, Jabbar*) 

fxo = 2.411046012096894, 
^1 = -0.5162699676752518, 
^2 = 0.15512703899481983, 
^(n = 2.1) = 2.360970285719317, 

(* 0.1153592 * 4.523926^ = 2.3609305957478286, Horedt*) 
fx{n = 1.9) = 2.4642242792543674, 

(* 0.140123 *4.1936152 = 2.4642600755839923, Horedt*). 
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10.2.6. n^5/2 

n = 5/2, 

Xq = 5.3552754590108239, 

^^(Xo) = -0.076264913479991799518, 

yl{Xo) = 0.0284821627061845798710092, 

t/Q"(Xo) = -0.0159555729247840799988, 
1/1 (Xo) = 0.1855760827385399754972, 
l/i(Xo) = 0.01349086250870339952973971, 
yl{XQ) = -0.005038344941156511, 
1/2 (Xo) = -0.03834028384550369027000, 
y'^(Xo) = -0.003769483809799692567, 

Xi = 2.43330 87690 08517 98883, (53) 
X2 = 1.03335 16315 07183 66045, 

X{n = 2.6) = 5.608939852226747, (* 5.60938 27386, Jabbar*) 
X{n = 2.4) = 5.122278098425044, (* 5.122278, Jabbar*) 
liQ = 2.18719956551707901061, 
= -0.3869041118787603560, 
fX2 = 0.1081049328559264189, 
;x(n = 2.6) = 2.14959 02036 57762, 
(* 0.06831578 * 5.609382^ = 2.14956 73869, Horedt *) 
f^{n = 2.4) = 2.2269710260 33514, 
(* 0.08489109* 5. 121870^ = 2.2269948490, Horedt*). 



10.2.7. n^3 



n — 3, 

Xq = 6.89684861937696, 
y'^^{Xo) = -0.04242975760445296, 
y'^{Xo) = 0.012304100016127638, 
y'^{Xo) = -0.005352053102589007, 
yi{Xo) = 0.16547682943945483, 
y[{Xo) = 0.006167361415567608, 
y'^{Xo) = -0.0017884578177459688, 
l/2(Xo) = -0.030535151278657177, 
y'^{Xo) = -0.0017139698191338721, 
Xi = 3.9000182603467954, 
X2 = 2.0525978144309858, 

X{n = 3.1) = 7.3073764235, (* 7.30848 42924 Jabbar*) 

X{n = 2.9) = ",6.52737277148659, (* 6.52637 41261Jabbar 

fio = 2.0182359509662287, 

/ii = ", -0.29335992553948076, 

^2 = 0.08152758134658877, 

^{n = 3.1) = ",1.9897152342 25746, 

(* 0.03725063 * 7.308484^ = 1.98970 28553, Horedt *) 

f^{n = 2.9) = 2.0483872193336428, 

(* 0.04809180 * 6.526374^ = 2.04840 08528, Horedt *). 



-28- 



10.2.8. n^7/2 

n = 7/2, 

= 9.53580534424479, 
yg(Xo) = -0.020790983939331436621, 
yl{Xo) = 0.00436061416708334168, 

^^"(Xo) = -0.00137186551413256, 
yi{Xo) = 0.148161582413238524348, 
y[{Xo) = 0.0024155338182516865, 
y'l{Xo) = -0.00050662397795474, 
y^{Xo) = -0.0249346780046498475, 
y'^{Xo) = -0.0007446596956576678, 

Xi = 7.12624197323115573338729, (55) 
X2 = 4.9541783719117226, 

X(n = 3.6) = 10.297971325287023, (* 10.30159 12890, Jabbar*) 

X{n = 3.4) = 8.872722930640792, (* 8.8695802788, Jabbar*) 

fio = 1.890557093443116455, 

/ii = -0.219648315244396768, 

H2 = 0.067713085342, 

pi(n = 3.6) = 1.86926 9392772097, 

(* 0.01761417* 10.301592 = 1.86926 42743, Horedt*) 

fx{n = 3.4) = 1.91319 90558 20976, 

(* 0.02431955 * 8.869580^ = 1.91320 56075, Horedt *) 
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10.2.9. n^4 



n = 4, 

Xo = 14.971546348838093, 
y'f^{Xo) = -0.008018078806403244, 
y'^{Xo) = 0.0010711089715893653, 

yQ{Xo) = -0.0002146289260907, 
1/1 (Xo) = 0.13250008777227387, 
y[{Xo) = 0.0006870175622715013, 
y'^{Xo) = -0.00009177643327735738, 
1/2 (Xo) = -0.020971140923665713, 
y'^lXo) = -0.0002923287696798706, 

Xi = 16.525166560655304, (56) 
X2 = 17.040461450049975, 

X(n = 4.1) = 16.794467619, (* 16.8137738705, Jabbar*) 
X{n = 3.9) = 13.48943 43072, (* 13.47384 13948, Jabbar*) 
fxo = 1.7972299144392503, 
fii = -0.1539930629857884, 
= 0.06552467522560379, 
f^{n = 4.1) = 1.78248 58548 92927, 
(* 0.006305183 * 16.81378^ = 1.78249 53973, Horedt*) 
f^{n = 3.9) = 1.8132844674 90085, 
(* 0.009988109 * 13.47384^ = 1.8132848994, Horedt*). 
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10.2.10. n^9/2 

n = 9/2, 

= 31.836463244695223, 
y^{Xo) = -0.0017145489124289180377, 
yl{Xo) = 0.0001077097602991190, 

yQ{Xo) = -0.000010149660105577, 
yi{Xo) = 0.1171972229781507365218, 
y[{Xo) = 0.000078867143804445547, 
y'^{Xo) = -4.95451666212243895495976 • 10"^ 
l/2(Xo) = -0.01868541777734960517998, 
y'^{Xo) = -0.00008663458543571578, 

Xi = 68.35455210905772991040946, (57) 
X2 = 139.00687095347173451822, 

X{n = 4.6) = 40.06198 71651, (* 40.41322 96343, Jabbar *) 
X{n = 4.4) = 26.3910767433, (* 26.15891 96485, Jabbar*) 
/io = 1.7377988676660324, 

= -0.0799366131848989555, 
H2 = 0.087809384368989, 
f^{n = 4.6) = 1.7306833001912325, 
(* 0.001059709 * 40.413242 = 1.73074 84954, Horedt *), 
pi(n = 4.4) = 1.7466706228282123, 
(* 0.002552517 * 26.15892^ = 1.74665 95493, Horedt *). 
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ll. Summciry 

In this paper we present the exact numerical solutions for the internal structure and 
global parameters of polytropic models of arbitrary index < n < 5. The perturbation 
method used here is not rigorously founded by means of the theory of differential equations 
and there are some problems about application of the method in the interval of argument 
where the perturbation function is of the same order or even larger than the initial non- 
perturbed function. The problem of justification is here similar to the problem of rotationally 
distorted polytropes which was already discussed in the astrophysical literature (see e.g. 
Chandrasekhar & Lebovitz (1962)). Still validity of any method in applicational sciences 
(as astrophysics is such relative to mathematics) may be checked with numerical calculations 
and (astro)-physical "common sense", and we showed in this paper that the perturbation 
method is applicable to the problem of the structure and the global parameters of the 
polytropic models. 
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